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ABSTRACT 

We explore the poloidal structure of two-dimensional magnetohydrodynamic (MHD) winds in relation to their 
potential association with the X-ray warm absorbers (WAs) and the highly ionized ultra-fast outflows (UFOs) in 
active galactic nuclei (AGNs), in a single unifying approach. We present the density n(r, 0), ionization parameter 
A u 0), and velocity structure v(r, 6) of such ionized winds for typical values of their fluid-to-magnetic flux ratio, F, 
and specific angular momentum, H, for which wind solutions become super-Alfvenic. We explore the geometrical 
shape of winds for different values of these parameters and delineate the values that produce the widest and narrowest 
opening angles of these winds, quantities necessary in the determination of the statistics of AGN obscuration. We 
find that winds with smaller H show a poloidal geometry of narrower opening angles with their Alfven surface 
at lower inclination angles and therefore they produce the highest line of sight (LoS) velocities for observers at 
higher latitudes with the respect to the disk plane. We further note a physical and spatial correlation between the 
X-ray WAs and UFOs that form along the same LoS to the observer but at different radii, r, and distinct values of 
n, £, and v consistent with the latest spectroscopic data of radio-quiet Seyfert galaxies. We also show that, at least 
in the case of 3C 111, the winds’ pressure is sufficient to contain the relativistic plasma responsible for its radio 
emission. Stratified MHD disk winds could therefore serve as a unique means to understand and unify the diverse 
AGN outflows. 

Key words: accretion, accretion disks - black hole physics - galaxies: active - quasars: absorption lines - X-rays: 
galaxies 
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1. INTRODUCTION 

Outflows are a common occurrence in accretion-powered 
objects. Their presence is affirmed by blueshifted absorption 
features in their optical, UV and X-ray spectra at frequencies 
of well defined transitions. Their velocities span a range of 
10 3 , roughly between 100 and 100,000 km s _1 , while the broad 
(~ 10,000 km s _ 1 ) absorption troughs of UV transitions are 
the defining characteristic of the class of broad absorption line 
quasars (BAL QSOs). The launch of Hubble, ASCA, Chandra, 
XMM-Newton, and Suzaku and their superior sensitivity, res- 
olution, and energy bandwidth showed that ~50% of active 
galactic nuclei (AGNs) exhibit signatures of such outflows in 
their optical, UV and/or X-ray spectra. Of these, the high res- 
olution long X-ray observations of Chandra and XMM-Newton 
are of particular interest because they discovered transitions 
that span a range of ~10 5 in ionization parameter, 7 /, (such as 
Fe i through Fe xxvi among others) and therefore sample a very 
broad range of conditions for the column density and velocity of 
plasma along the observers’ line of sight (LoS). The fact that the 
plethora of these transitions spanning 5 decades in / is squeezed 
within roughly 15 years of frequency underscores the utility of 
X-ray spectroscopy. 


7 f = L/(nr 2 ) where L is an ionizing luminosity (usually defined between 1 
and 1000 Ryd), n is the plasma number density, and r is distance from the 
ionizing source. 


The first signature of X-ray absorbing plasma in AGN was 
that in th e. Einstein spectrum of the QSOMR 2251-178 (Halpern 
1984), attributed to “warm” plasma of temperature T ~ 10 6 K 
(rather than “cool” T — 10 4 K clouds), thus coining the term 
warm absorber (WA) for absorption features in the <1 keV 
band. The ubiquity of WAs was eventually established by ASCA 
which discovered absorption features of typical column density 
of Nft ~ 10 2 °-10 22 cm -2 and £ ~ 10 _1 -10 4 erg cm s _I at 
moderate outflow velocities (v < 3000 km s -1 ) in the spectra 
of ~50% of Seyfert 1 galaxies (Reynolds 1997; George et al. 
1998). Lollowing these discoveries, extensive spectroscopic 
observations with Chandra and XMM-Newton have been made 
to study in detail the physical conditions and spatial origin of 
WAs (e.g., Blustin et al. 2005; McKernan et al. 2007; Reeves 
et al. 2009a; Turner & Miller 2009; Torresi et al. 2010, 2012). 
These are of interest because the presence of ionized absorbing 
gas of a very wide range of / along the observers’ LoS offers 
the opportunity to map the spatial distribution of this gas: one 
can easily see that the knowledge of / (from the presence of 
certain ionic species) and the measurement of their hydrogen 
equivalent column Ah [< n(r)r] in the observed spectrum can 
provide a measure of the absorbing gas’ density dependence 
on the distance from the continuum source along the observer’s 
LoS. Such information is extremely valuable in assessing the 
global properties of the winds which presumably give rise to the 
observed absorbers. 

A quantitative, more systematic formal way of implementing 
the above procedure is to construct the absorption measure 
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distribution (AMD), namely the hydrogen-equivalent column of 
specific ions per decade of f ; i.e., AMD = clNyi/dlog^ (e.g., 
Holczer et al. 2007; Behar 2009; Detmers et al. 201 1; Holczer & 
Behar 2012, among others). These authors, instead of adding gas 
components at different values of f and Ah until a satisfactory 
statistical significance x 2 is obtained, assumed a continuous 
dependence of binned Ah on £ (i.e., c/An/t/logl oc f 1 ) and 
provided a global fit to the entire set of ionic transitions with 
s as the sole parameter. The continuous dependence of Ah on 
f also implies (assuming a smooth spatial gas distribution) a 
continuous dependence of the ionized medium density with r, its 
distance from the AGNs, of the form n(r) oc r -( 2s+1 )/(s+D X-ray 
analysis of the spectra from a number of Seyfert galaxies showed 
that Ah has only a weak dependence on £ (s — 0) (Behar 2009), 
implying a wind density profile of n(r) oc r ~ 1 (the largest value 
of s found was .v ~ 0.3 implying n(r) oc r _L25 ). Furthermore, 
the presence of UV absorption features in the spectra of AGNs 
that exhibit X-ray WAs (e.g., Crenshaw et al. 2003) suggests a 
physical link between these components, which at present is not 
completely understood. 

In parallel with the high-resolution observations discussed 
above, a number of lower resolution CCD X-ray observations 
of radio-quiet Seyferts and BAL QSOs have discovered ab- 
sorption features identified as highly ionized iron (primarily 
Fe xxv/Fe xxvi at energies ~7-8 keV in the source frame) of 
column densities (Ah ~ 10 23 -10 24 cm~ 2 ), generally higher than 
those of their moderately ionized ions at higher (blueshifted) 
velocities ( v/c ~ 0. 1-0.7 where c is the speed of light), called 
ultra-fast outflows (UFOs; Tombesi et al. 2010a). These are 
seen across different AGN populations such as Seyferts, BAL 
QSOs and non-BAL QSOs (Chartas et al. 2003, 2007; Reeves 
et al. 2009b; Pounds & Page 2006; Tombesi et al. 2010a, 201 la, 
2012a). It is worth noting that similar features (UFOs) have 
also been reported in radio-loud Seyferts such as 3C 111 (e.g., 
Tombesi et al. 2010b, 2011b), indicating that these high outflow 
velocity transitions represent a generic feature of AGN rather 
than one associated with a specific AGN subclass. Furthermore, 
in the case of a BAL QSO APM 08279+5255, Chartas et al. 
(2009) noted a likely correlation between X-ray photon index 
r and the measured outflow velocity of Fe xxv, implying ad- 
ditional underlying physics pertaining to the velocities of these 
absorbers. 

These observations of high ionization (logf > 4), high 
velocity (v 2; 0.1c) outflowing gas (UFOs) are of importance 
because they challenge the conventional radiative acceleration 
scenario of AGN outflows (e.g., Proga et al. 2000) which 
demands low ionization for this gas so that line radiation 
pressure driving is efficient; as such, they lend support to 
magnetic driving for these flows. The recent discovery of the 
UFOs in radio-loud Seyferts exhibiting jets (Tombesi et al. 
2010b, 2011b) also seems to challenge a simplistic wind-jet 
dichotomy that associates “winds” with radio-quiet AGNs and 
“jets” with the radio-loud AGNs, perhaps implying a multi- 
component gas in the same outflowing plasma across these 
broad AGN classes. While few simultaneous X-ray-radio studies 
having been conducted so far in radio-loud AGNs, Tombesi et al. 
(2012c, hereafter T 12c) performed such a combined study of the 
X-ray UFOs and the Very Long Baseline Array (VLBA) radio 
jet of the radio galaxy 3C 111; their simultaneous presence in 
this object has indicated that these components can actually 
coexist in AGNs. Also, with respect to the coexistence of WAs 
and UFOs in radio-quiet Seyferts (e.g., Tombesi et al. 2010a, 
2011a, 2012a, 2012b. hereafter T 12b), have analyzed X-ray data 


from 35 radio-quiet Seyferts to obtain a number of correlations 
between f , Ah, and v of these components that argue in support 
of a single physical entity underlying the nature of WAs and 
UFOs. 

Similar absorption features have been observed and studied 
intensively in binary systems including transient black hole 
(BH) candidates and neutron star low-mass X-ray binaries. 
These typically exhibit highly ionized plasma (e.g.. Si, S, Mg, 
Ar, and Fe) of log § ~ 3-4 and A H ~ 1 0 2 1 — 1 0 22 cm -2 at typical 
outflow velocities ranging from v ~ 300 km s~ 1 to 2000 km s _ 1 , 
i.e., properties similar to those of AGNs but with lower velocities 
and more highly ionized species. Recently, however, King et al. 
(2012) found an extremely fast (v/c ~ 0.03-0.05) ionized 
outflow based on the Fe xxv K-shell signature from Suzaku/XIS 
observations of a transient X-ray source, IGR J17091-3624. 
This is currently the fastest X-ray wind known in galactic 
binaries. This clearly illustrates a broad range of X-ray outflow 
kinematics possible around stellar-size compact objects. Some 
of the BH binaries, such as GRS 1915+105 (e.g., Neilsen & 
Lee 2009; Neilsen et al. 2012) and H 1743-322 (e.g., Miller 
et al. 2006; Blum et al. 2010), exhibit relativistic radio jets 
akin to those of radio-loud Seyferts. In particular, these two 
components of outflows. X-ray absorbers and radio jets, appear 
to be anti-correlated: X-ray winds seem to be suppressed in 
the low/hard state of these objects which favor the presence of 
radio jets (e.g., Fender et al. 2004). Also, VLBA monitoring of 
jet emission combined with simultaneous RXTE observations of 
low-frequency quasi-periodic oscillations seems to imply that 
the jet ejection is anti-correlated with the disk activity (see, e.g., 
Miller- Jones et al. 2012 for utilizing QPOs as a proxy of disk 
activity in H 1743-322). These observations are suggestive of 
a mutual interplay between the disk (where X-ray outflows are 
also launched) and jet production. 

In an effort to place the ensemble of these facts into the 
broader unifying context of accretion flows onto compact 
objects and their associated winds, Fukumura et al. (2010a, 
hereafter FKCBlOa) have considered the photoionization of 
magnetohydrodynamic (MHD) accretion-disk winds. Motivated 
by the observed broad range of the WA and UFO ionization 
parameter and velocity range, they opted for self-similar models 
of MHD winds, as they naturally cover many decades in radius. 
Since their enunciation by Blandford & Payne (1982, hereafter 
BP82), these winds have been the subject of many studies, 
both semi-analytical (BP82; Contopoulos & Lovelace 1994, 
hereafter CL94; Konigl & Kartje 1994, hereafter KK94; Everett 
2005; Ferreira 1997) and fully numerical (Proga 2003; Pudritz 
et al. 2006; Fendt 2006; Ohsuga et al. 2009; Porth & Fendt 
2010; Murphy et al. 2010). A crucial hint for the choice of 
a specific model has been the value of the parameter s(— 0) 
in the AMD relation (the dependence of Ah on §), since this 
determines the run of the wind density with distance. Clearly, 
observations favor density profiles close to n(r) oc 1/r, a 
value significantly different from that of the winds of BP82 
(which imply s ~ 3 /2) and at significant odds with radiation- 
driven outlflows (which, for / approaching its asymptotic value 
foo, imply .s 1 and hence n(r) oc r~ 2 ). The important 
features of these density profiles are their almost constant 
(or slowly decreasing) column density per decade of radius 
and their linear decrease of £ with distance r, a fact that 
determines that high ionization ions have higher velocities than 
the lower ionization ones, in agreement with observations. The 
models of FKCB 1 0a reproduce these features yielding velocities 
v ~ 300 km s’ 1 for the moderately ionized species of Fexvn 
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and v ~ 3000 km s _I for the highly ionized ions of Fexxvi. 
This work was extended by Fukumura et al. (2010b, hereafter 
FKCB 10b), which demonstrated the model’s broad applicability 
by successfully modeling the combined X-ray UFOs and UV 
broad absorption features of the BAL QSO APM 08279+5255 
(Chartas et al. 2009); this was achieved by adjusting only the 
relative normalization of the UV and X-ray fluxes to those 
appropriate for BAL QSOs. 

To date, our publications on the subject (FKCB 10a, 
FKCB 10b), which presented the first attempt in this specific 
direction, have employed a single set of wind parameters im- 
plying an identical magnetic field geometry in both cases. The 
goal of the present work is to remedy this deficiency by an explo- 
ration of the wind geometry and ionization through a parameter 
search, keeping in mind the X-ray and radio AGN observations 
(T 12 b; T 12 c), with the underlying theoretical framework still 
based on the work of CL94. In Section 2 we summarize the 
basics of MHD accretion-disk wind models that have been pre- 
viously applied for comparison with the observed X-ray WAs 
(FKCB 10a) and UFOs (FKCB 10b). In Section 3 we present our 
results and demonstrate a relevance of the outflows in a num- 
ber of situations where those winds could provide a substantial 
impact on the surrounding environment, and we discuss the im- 
plications of our results in Section 4 with a comparison with 
BP82 winds in the Appendix. 


2. CHARACTERISTICS OF MAGNETIZED 
ACCRETION-DISK WINDS 


As is well known, the structure of axisymmetric MHD flows 
is determined by the magnetic flux function 'P(r, 0 ) (see, e.g., 
CL94, Ferreira 1997; Vlahakis & Konigl 2003; Fendt 2006; 
Porth & Fendt 2010). As discussed earlier, the broad range of 
the observed £ and v of absorption features suggests a power- 
law form of all physical quantities with the radial coordinate r; 
as a result, we assume for T^r, 0 ) the following separation of 
variables: v P(r, 0) in the power-law form of spherical radius r as 
v P(r, 0) = (r / r o ) qx V{0yV 0 , where q is a self-similar index and 
V P 0 is the magnetic flux normalization at the innermost radius 
r a at which our scaling is assumed to be valid, and M'( 0 ), the 
angular dependence of this function which has to be determined 
numerically. We provisionally assume r — r a = 3Rs ( Rs is the 
Schwarzschild radius), the value of the innermost stable circular 
orbit of the Schwarzschild geometry. We make this choice based 
on the high velocities (u ~ 0.5c) observed in certain objects, 
with the understanding that those of our results pertaining to 
these innermost radii may not be accurate. However, we do not 
expect this assumption to affect the wind structure at the much 
larger scales covered by the observations. 

With the above choice for V F the magnetic field strength, 
velocity, number density, and total pressure take the following 
form 


\B(r,0)\ = B 0 

|t)(r, 0 )| = v 0 ( 

1 { B„ 
n(r , 0 ) = — — 

m p \v 0 

p(r , e) = KBl I 


£ 

2 


r 

r a 


B(9), 

2 

0 ) 

m, 

( 2 ) 

2q-3 


) m. 

(3) 

■4 


m r , 

(4) 


where 


m 


Fq B p (9) 
4tt v p (9) ' 


(5) 


m p is proton mass, and K is the dimensionless polytropic con- 
stant related to plasma entropy. All the dimensionless angular- 
dependent functions denoted by tildes must be numerically ob- 
tained from the conservation equations as the solution to the 
Grad-Shafranov equation with initial values on the disk at 
(r, 6) = (r a , 90°) to which all the quantities are normalized 
as denoted by the subscript “o.” The magnetized wind is also 
characterized by its plasma /3, defined as the ratio of thermal to 
magnetic field pressures 


Mr, 9) = m = 


P(r , 0)/K 
B(r, 0 ) 2 


n(9) r 

lm 2 ’ 


( 6 ) 


where T = 5/3 is the polytropic index of the wind plasma. We 
note that its LoS radial dependence drops out in our framework 
since B(r) oc r q ~ 2 and p(r) oc r 2q ~ 4 . 

Under steady-state axisymmetric conditions there are five 
conserved quantities along a streamline of given v P(r, 0 ); i.e., 
the particle flux to magnetic flux ratio F, the angular velocity of 
field lines £2, the total angular momentum of the plasma H, the 
total energy (Bernoulli function) J, and the entropy S (see CL94 
for details). F and H are defined as 

IvJ 

FC¥) = 4nnm D F- E F, (7) 


H(¥)= -F(V)rv t +rB t , ( 8 ) 

where v p and B p are, respectively, the poloidal velocity and 
magnetic field while v$ and B,i, are, respectively, the azimuthal 
velocity component and magnetic field strength. 

We assume the disk threaded by the magnetic field to 
be infinitely thin, so we define boundary conditions at the 
point ( r = r a ,9 = 90°) where v = v„ ~ with only 
a small component v Zt0 perpendicular to the disk surface 
[iy(90°), Vo(90 ) v r i,]. For a given set of F a and //„, a trans- 
Alfven wind solution can be obtained when the initial magnetic 
field line orientation 9 0 is met with the regularity condition . 8 
With these constraints, and J B are dependent quantities 
constrained by F a and H a (see CL94 and FKCB 10a); i.e., 

a, = 1 - v e (90°)(F o + H 0 ), (9) 


Jo = Vg(90°)[l + if{ 90°)2]/2 - 1/2 - £ 2 C , ( 10 ) 

where i//-(90°) characterizes the angle the poloidal field/flow 
lines make with the disk at their footpoints (readers can 
easily verify that this angle is equal to tan~' {< 7 /t/t( 90°)} with 
vcp( 90°) = 1 and vg( 90°) = 0.01 at the launching radius 
r 0 (see CL94 for details). We thus find that, for a given 
q, the wind structure is primarily governed by F a and H a ; 
the goal of the present study is to determine the structure 
and observable properties (i.e., column density, LoS velocity 
as a function of angle) of these winds and their relation 


8 While physically valid winds must pass through a slow magnetosonic point, 
we assume that the wind quickly becomes super-slow magnetosonic due to its 
efficient acceleration (see CL94 and Ferreira & Pelletier 1995 in relation to the 
wind conditions). 


3 


The Astrophysical Journal, 780:120 (11 pp), 2014 January 10 


Fukumura et al. 


to observations. Specifically, the value of F 0 is generally 
responsible for determining the wind kinematics while H a plays 
an important role in shaping its poloidal structure. Note that one 
must adjust either F a or H 0 for a given initial poloidal angle of 
the magnetic field in order for winds to become super- Alfvenic 
(see the Appendix for a discussion on the viability for our choice 
of parameter sets). Although the other model parameters (see 
CL94 and FKCBlOa) are also coupled to and affect the wind 
properties, we find them to be fundamentally less significant 
compared to F 0 and H„. 

The magnetic and velocity fields can be decomposed explic- 
itly in poloidal and toroidal components as B(r. 0) = B p + B,/, 
and v(r, 0) = r sin QQ.<j) + m p F /(A tzti)B, while the density at 
(. r a , 90°) can be explicitly expressed in terms of the dimension- 
less mass-accretion rate m (see FKCB 10a) as n 0 = n(r 0 . 90°) = 
rtw>h/(2(7Trs), where ijw is the ratio of the total mass flux rate 
launched in the wind to m (while the rest is accreted) and ay is 
the Thomson cross section. Note also that in our discussion on 
density normalization that m always refers to the local mass flux 
at the innermost flow radius at r ~ r a because the wind mass 
flux generally has a radial dependence. One can also calculate 
force elements in the wind (e.g., gas pressure force and mag- 
netic forces) as V p w j n d = V p gas — (J x B)/c (where J is the 
current density of the wind) which can be a measure of thermal 
equilibrium with ambient gas. 

We now briefly discuss the issue of critical points in our winds. 
These are discussed in BP82, CL94, and also in Vlahakis et al. 
(2000) and Ferreira & Casse (2004). A complete description 
would require that these solutions pass the fast and modified fast 
(i.e., the critical point in the direction perpendicular to that of the 
imposed symmetry) points as described in BP82 (and achieved 
in Vlahakis et al. 2000; Ferreira & Casse 2004). However, these 
points lie at sufficiently high latitudes where the wind densities 
are very low and the gas ionization complete, so that they are 
effectively irrelevant to the absorption features observations. As 
a result, we limit our wind description to crossing the Alfven 
and the slow-magnetosonic points. While crossing all critical 
points is an issue of principle with steady-state winds (e.g., 
BP82 and CL94 self-similar winds recollimate on the axis and 
even turn back), these points are in fact crossed without much 
trouble in nature; i.e., in time-dependent simulations (Pudritz 
et al. 2006; Fendt 2006; Porth & Fendt 2010). Once the fluid 
has achieved escape velocity, it reaches infinity (albeit after a 
potential recollimation along the symmetry axis). 

In the following calculations discussed in Section 3, therefore, 
we chose to stop the wind solutions past the Alfven point 
arbitrarily to eliminate a physically less reliable portion of it. 
However, since our primary goal has only been to demonstrate 
the model’s global viability with the observed X-ray absorbers 
and not to construct a comprehensive model of ultrarelativistic 
winds/jets, this approximation is justified up to near the Alfven 
point. Most notably, we find that this manipulation makes little 
impact on our argument of WAs and UFOs in the context 
of MHD-driven outflows presented in this paper since most 
absorbers still lie within such an artificial cut-off. 

3. RESULTS 


Table 1 

Characteristics of Baseline MHD Wind Models 


Model 

Mr) 

F„ 

H 0 

0o 


Jo 

e A 

v p a / v 0 

$open,4 

A 

r~ l 

0.065 

-1.7 

1.276 

1.016 

-1.51 

38° 

4 

7° 

B 

r- 1 

0.03 

-2.455 

1.872 

1.024 

-1.52 

53° 

7 

28° 

C 

,.-1.5 

0.05 

-3.152 

5.7 

1.032 

-1.53 

<1 

ON 

o 

6 

47° 

D 

,.-1.5 

0.1 

-1.68 

2.105 

1.015 

-1.515 

59° 

3 

6° 


they are best known and most widely cited. Table 1 shows the 
functional form of «(r), the values of F a , ll 0 , il„, and J a , the 
angle of the Alfven surface 0 a, the wind asymptotic (i.e., at 
z/r c = 10 4 ) poloidal velocity normalized to v 0 , v p a/v 0 , and 
the wind opening angle 0 op en ,4 at z/r 0 = 10 4 for the wind 
element launched at the fiducial radius r 0 . The corresponding 
poloidal wind structure of all four cases is shown in Figure 1 
in linear scale in the computational domain 1— 10 4 r o . The 
figure shows the poloidal density distribution normalized to 
its maximum value n a = B 2 /(m p v 2 ) (in fact it shows the 
value of log[m p n/(B 0 /v 0 ) 2 ]) along with its contours (solid thin 
lines), the wind streamlines r(6) (solid thick lines) equally 
spaced by A r/r a = 10 3 , the poloidal velocity vectors v p /v 0 
(gray arrows), and the velocity contours (dashed lines). A white 
solid line denotes the Alfven surface where the plasma poloidal 
speed v p reaches the local Alfven speed v a = B p / \J 4nm p n. 
As mentioned earlier in Section 2, we decided to artificially 
truncate the super-Alfven wind solution at a critical latitude 
0 C = cot - '(100 cot 6a) because of uncertainty in self-similarity 
as well as trans-fast nature of the solutions where 0 a denotes 
the latitude corresponding to the Alfven point. 

The values of variables in these figures in physical units 
can be obtained from Equations (l)-(6) with an appropriate 
normalization for the plasma density and velocity. While the 
velocity normalization can be obtained independently of the 
mass of the accreting object M, once a choice of the fiducial 
radius r a (close to the Schwarzschild radius Rs) is made (where 
v 0 — c (r 0 / R S F^ 2 ), the normalization of the density (and the 
ensuing ones of the pressure and magnetic field) requires a 
value for the BH mass M and the dimensionless accretion rate 
m (normalized to its Eddington value). The normalization of 
the density can be obtained by considering that for accretion 
or outflow at the local escape velocity and at the Eddington 
rate, the Thomson optical depth of the flow is equal to one 
or greater at r a ~ Rs', this implies that n a — iii/a T R s or 
n a ~ 5 x I O 10 (m/ Ms) cm -3 . With these considerations and 
assuming magnetic field in equipartition with the kinetic energy 
density, i.e., nm p v 2 /2 = B 2 /%n, in a flow with m ~ 1 and 
M — 10 8 M q we obtain the values 



Employing the framework discussed above, we have studied 
the structure and geometry of winds with two distinct density 
profiles, namely, winds with n(r) oc r -1 (i.e., q ~ 1; CL94, 
KK94), because they are favored by the X-ray absorber data, 
and winds with n(r ) oc r -3 / 2 (i.e., q ~ 3/4; BP82), because 


p(r,d) — 2.1 x 10 4 




/ m 

\n(90") 


r 

dyne cm -2 , 



(13) 


4 


The Astrophysical Journal, 780:120 (1 lpp), 2014 January 10 


Fukumura et al. 



I 


C/5 

c 

<L> 

O 


(r/ro)sin# 


I 


10000 



I 


C/5 

a 

o 

Q 


0 2000 4000 6000 8000 10000 

(r/ro)sin0 


I 


10000 


c 

'O 



I 


I & 


c 

« 

Q 


0 2000 4000 6000 8000 10000 

(r/ro)sin0 |_ I2 



I 


□ 

Q 


(r/ro)sin# 


I 


Figure 1. Dimensionless quantities of the MHD winds in the poloidal plane (r sin0, r cosQ) showing density distribution \og[m p n(r, 0)/(B o /v o ) z \ along with its 
contours (solid) superimposed by wind velocity v p (r, 0)/v o (arrows), its contours (dashed), and streamline geometry r(6) (thick solid) for the models (A)-(D). A 
white line denotes the Alfven surface. Distance is normalized to r Q around a BH. The density is normalized to its maximum value. The wind parameters in each model 
are listed in Table 1 . 

(A color version of this figure is available in the online journal.) 


with the angular parts of h, p, and B set equal to one on the 
disk surface ( 9 — 90°) and their 0-dependence produced by 
the solution of the MHD equations. The normalization of the 
quantities given above, appropriate for M = 10 8 M 0 , scales like 
1/M, i.e., inversely with the BH mass, while the plasma /I is 
independent of it. 

The properties of the models of Figure 1 are summarized in 
Table 1 (see the Appendix for a brief comparison with BP82 
wind parameters). All these models have approximately the 
same values of and J a \ models (A) and (B) have density 
profiles of n(r) oc r 1 (CL94; KK94), while models (C) and 
(D) have n(r) oc r~ 3 ^ 2 (BP82). The larger value of H„ in 
model (B) results in less collimation, i.e., larger asymptotic 
opening angle (0 O pen,4 ~ 28°), compared with that of model (A) 
($open .4 ~ 7°) which has a smaller H„ value. For the same reason, 
the Alfven surface of (A) lies at smaller value of 9 (9a ~ 38°) 
than of (B) (9a ~ 53°). Finally, the combination of larger H a 
and smaller mass loading leads to larger asymptotic value in 
(B) (u p , 4 /u 0 ~ 7) compared to that of (A) (v v a/v 0 ~ 4). The 
models with n(r) oc r -3 / 2 (C,D) have relative dependencies 
similar to those of (A, B) and are summarized in Table 1. 
Finally, for models with similar values of J 0 , and H„ 
but different density profiles, namely (A) and (D), while the 
asymptotic opening angles at z = 1 0 4 r o are similar, the Alfven 


angle is smaller in (A) than (D) presumably as the result of the 
different poloidal current distribution of these two solutions via 
a coupling effect between F a and H a . 

These results are summarized in Figure 2(a) which shows 
the outflow opening angle 0 open along the innermost streamline, 
originating at r = r a (these curves correspond to the innermost 
streamlines in Figure 1) as a function of axial distance z/r a 
for models (A)-(D) with dots denoting the Alfven points. The 
effects of the higher H a in models (B) and (C) are apparent there. 

Our interest in the driving mechanism focuses primarily on 
the J x B force while not discussing the centrifugal force in 
detail. In the wind model considered here, toroidal wind velocity 
scales as v,/, oc r 1/ 2 along a given LoS. Therefore, centrifugal 
force term can be written as ~pv^/r oc r~ 3 along the LoS 
and magnetic force term also can be written as VB 2 oc r~ 3 for 
n oc r _1 . The relative strength of the centrifugal term can be 
dominant at smaller radii at least within the Alfven surface. This 
behavior is indeed confirmed in Figure 2(b), which compares 
the force components tangent to the poloidal field line (i.e., 
accelerating components) following Lovelace & Romanova 
(2013). The relative strength of each force component and its 
profile is consistent with those in other works (e.g.. Figure 9 
in Porth & Fendt 2010) in that the magnetic part can be the 
dominant component even at smaller radii. Lii et al. (2012) also 
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Figure 2. (a) Disk- wind opening angle 0 O pen as a function of axial distance z (up to 6 C ) along the innermost streamline (with 4* = const) of winds launched at r = r a 
for models (A)-(D). (b) Force components tangent to the poloidal field line for model (A) showing J x B term Fm (dashed), centrifugal term F c (solid) and gas 
pressure term F p (dotted). The wind parameters are listed in Table 1. 


argues in their 2.5D MHD jet simulations that the centrifugal 
force is cancelled by gravity and thus the jet is driven by a purely 
magnetic force. This mechanism is similar to the inner disk 
wind model discussed in Lovelace et al. (1991) and observed in 
simulations of conical winds in Romanova et al. (2009). 

Focusing more on the global properties of magnetized wind, 
Figure 3 shows kinematic profiles v of the wind and the magnetic 
field B as well as the wind density profile n as a function of 
axial distance z/r a for model (A) in the upper panels and (B) 
in the lower panels. Due to a geometrical difference between 
the two models, the radial wind velocity still exceeds the axial 
one in model (C). While velocity profiles v are independent 
of the models since we assume Keplerian boundary condition 
at the disk surface, the magnetic field profiles B and density 
n are indeed dropping faster (i.e., oc z" 3 / 2 ) in model (C) as 
expected. It is seen in both models that the wind undergoes a 
rapid acceleration phase from the base through the Alfven point 
(i.e., the intersection between v z and va ), after which the wind 
asymptotically approaches a coasting speed at v v a/v„ ~ 4 in 
model (A). The predominant initial orbital motion ( v,j, ) of the 
wind at the time of launch is found to be efficiently converted 
into axial motion (v-). These profiles follow |u| oc z" 1 / 2 , |Z?| oc 
z~ l and n oc z 1 / 2 at large distances in consistence with the 
chosen 1/r self-similarity in model (A). 

As discussed in Section 1, the observations and physics of 
WAs mainly provide their ionization parameter £ and hydro- 
gen equivalent column density Ah, and, to the extent that 
can be measured, also the plasma velocity along the ob- 
server’s LoS. Therefore, in order to relate the above mod- 
els to observations, we plot in Figure 4 (in log-log space) 
the structure of ionization parameter log £ [erg cm s" 1 ] for 
m = 0.1, along with the normalized density contours (dot- 
ted curves for 10" ', 10" 2 , 10" 3 , 10“ 4 , 10" 5 , 10" 6 , 10" 7 from 
bottom to top) and velocity contours (solid black curves for 
v/c = 0.5, 0.3, 0.1, 0.05, 0.03, 0.01 from innermost to outer- 
most) of the same winds shown in Figure 1 . Here, we only con- 
sider the geometrical dilution factor (see footnote 1), whereas 
in our photoionization calculations the opacity of the wind 
material was included (FKCBlOa, FKCBlOb). Because these 
quantities span a range of several decades, it is only reason- 
able to present the wind structure in logarithmic coordinates. 
Thus, color-coded log £ is shown in the range ~0 and ~10. To 
make contact with observations, we also draw the observer LoS 


of different inclinations (diagonal gray lines corresponding to 
0 — 10°, 30°, 45°, 60°, 80° from top to bottom). One can now 
see that low inclination lines of sight intercept mainly high ion- 
ization and low column plasma (the iso-density curves are more 
closely spaced at the lower inclination directions). It is also ap- 
parent that the wind is highly ionized near its innermost edge 
where the velocities are also higher but the column intercepted 
is lower. It is not clear whether these sections of the wind can be 
detected in actual observations. With increasing inclination the 
LoS samples lower values of f and velocity v and higher values 
of the column Ah- This model therefore naturally associates 
the observed WAs with objects of sufficiently high inclination 
angles. 

The two bottom panels of Figure 4 correspond to winds with 
n(r) oc r 3/2 (models C, D). These panels look quite similar to 
those of models (A) and (B) above them. The obvious difference 
is the more closely spaced iso-density contours due to the steeper 
density profile; as a result, most of the column along a given LoS 
is concentrated near log(r) = 0. For the cases of both density 
profiles, we see that the higher angular momenta solutions leave 
a large fraction of the LoS not intercepting any plasma, at least 
within the computational domain. This fact would be of interest 
for the statistics of absorbers. 

The column of these figures along a given LoS, one of the fun- 
damental observables of the AGN absorption features, depends, 
as discussed above, on the value of (the normalized) mass flux 
and, as is apparent in Figure 4, on the winds’ poloidal struc- 
ture. Clearly models (B) and (C) require observer inclinations 
greater than 45° to produce any significant absorption column. 
For model (A), a column of Ah — 10 23 cm" 2 and velocity 
~10, 000-15, 000 km s" 1 , consistent with the observations of 
UFOs, assuming mass flux m ~ 1 , demands an inclination angle 
6 ~ 30°-45°, based on the density function n(0) of FKCBlOa, 
while the observed velocity implies r/r a ~ (v/c)~ 2 — 10 3 . 
In model (A) we see that in this region of parameter space 
log £ ~ 3—4, consistent with the presence of an Fe xxv UFO. 
Clearly models (B) and (C) have very little column for incli- 
nations 9 & 30°. These figures also allow one to estimate the 
the velocity of a given ion (i.e., a given value of f ) at a specific 
LoS and as such they are of value in setting a particular set of 
observations within the framework of these winds. 

Besides the above connection between WAs and UFOs, 
we would further like to point out that the highly ionized 


6 


The Astrophysical Journal, 780:120 (1 lpp), 2014 January 10 


Fukumura et al. 



Figure 3. Fiducial wind velocity of radial v r , toroidal v axial v z components along with the poloidal Alfven speed va and slow magnetosonic speed normalized to 
v 0 (along the innermost streamline) and the corresponding magnetic field of radial B r , toroidal B $ and axial B z components normalized to B a as well as density profile 
n as a function of axial distance z/r c . We show this for model (A) in the upper panels and model (C) in the lower panels. The wind parameters are listed in Table 1. 


section of the wind near its axis provides thermodynamic 
characteristics consistent with radio jet properties. To this end, 
we show in Figure 5(a) (in linear space again), the innermost 
(r < 800r o ) structure of model (A): the contours of wind 
poloidal velocity v p /v D are given by the solid curves (for 
v p /c = 0.2, 0.1, 0.05, 0.03, 0.01 from innermost to outermost) 
while those of the wind pressure by the dashed ones (for 
Pueo = 10 _1 , 10 -2 , 10 — 3 , 1 0 4 dyne cm -2 from bottom to top) 
and are superimposed on the color coded ionization structure 
log |(r, 0), normalized to its maximum value. In this region, 
the pressure p(200r o , 40°) ~ 0.012 dyne cm -2 and velocity 
Vp/c > 0.01 are consistent with those obtained for the UFOs 
of the radio galaxy 3C 111 p > 0.001-0.01 dyne cm -2 
(T12c), which were also shown to be compatible with the 
pressure of radio jets in the same region. Because of their 
tighter collimation, we would therefore associate models (A) 
and (D) with the structure of radio-loud Seyferts. Finally, a 
choice between these two models will rely on the dependence 
of the value of Ah on the ionization parameter £ , assuming that 
the entire X-ray spectrum is produced by a point-like source, 
i.e., a source whose extent does not depend on the X-ray energy. 

4. SUMMARY AND DISCUSSION 

We have presented above the two-dimensional ionization 
structure of the MHD winds which our previous works we have 
associated with the WAs of AGNs. We have argued further 


that the same winds are also compatible with the more recently 
discovered UFOs and that in fact these features are different 
facets of the same underlying phenomenon, which apparently 
also occurs in radio galaxies, such as 3C 111. We have also 
argued that the radio/X-ray observations of this later object 
imply pressure equilibrium between its radio jet emitting plasma 
and the plasma associated with its UFOs, suggesting that the 
former occupies, in pressure equilibrium, the near-axis region 
of the self-similar UFO-WA flows of our models. 

The main results of the present work are encapsulated in 
Figures 1 and 4: they present the extreme cases of wind 
collimation geometry compatible with our assumptions for 
two different values of the wind radial density dependence. 
We found that while the density profile does affect the wind 
collimation, the latter is determined mainly by the value of their 
specific angular momentum H a , with higher values producing 
significantly less collimation than lower ones. The issue of jet/ 
MHD wind collimation is of interest in studies of AGN statistics, 
considering that the outer regions of these winds could play the 
role of the AGN unification “torus” (Konigl & Kartje 1994). 
Thus, higher collimation implies a larger fraction of sources 
with X-ray absorption features. A more collimated flow due to 
less ambient matter entrainment, is more likely to produce a well 
defined, large-scale narrow jet and is therefore of higher velocity. 
Similarly, less collimated flows are likely to entrain more 
matter leading to lower velocity jets. The point of the present 
work is to note that, while it is difficult to image the regions 
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Figure 4. Contour plot (in log-log space) of models (A)-(D) corresponding to Figure 1 showing LoS poloidal velocity (solid) of v/c = 0.5, 0.3, 0.1, 0.05, 0.03, 0.01 
(decreasing outward) and the density (dashed) log ft = — 1, —2, —3, —4, —5, —6, —7 (decreasing upward) superimposed by color-coded ionization parameter log£ 
which is normalized to their maximum value for m = 0.1. Five LoS angles are denoted by thick gray lines ( 6 = 10°, 30°, 45°, 60°, 80° from top to bottom). The wind 
parameters are listed in Table 1 . 

(A color version of this figure is available in the online journal.) 



Figure 5. (a) Contour plot of model (A) showing LoS poloidal velocity (from top to bottom in solid) of v/c = 0.2, 0.1, 0.05, 0.03, 0.01 and wind pressure (from top 
to bottom in dashed) of p\jfo = 10 _1 , 10 -2 , 10~ 3 , 10~ 4 dyne cm -2 superimposed on color-coded ionization parameter log £(r, 6) that is normalized to its maximum 
value, (b) A schematic description of stratified MHD-driven wind structure manifesting WAs and UFOs (not to scale) with column distribution (similar to Figure 8 in 
K12 and Figure 5 in T12b). The arrows indicate a mutual interaction between inner relativistic jets and the MHD winds at their interface. 

(A color version of this figure is available in the online journal.) 
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where these flows originate and collimate, one can obtain their 
physical conditions (column, velocity) from spectroscopic X- 
ray observations and also its collimation properties from the 
statistics of AGNs with a given amount of obscuration (Malizia 
et al. 2009; Terashima 2010) and their relation to the direction 
of a larger scale radio jet. A resolution of these issues, as well 
as those of systematics such as the obscuration dependence on 
luminosity (Tueller et al. 2008; Bourlon et al. 2011), require 
further observations and modeling (Proga 2003; Everett 2005; 
Kazanas et al. 2012). However, we believe that large-scale 
MHD winds will form the basis for an account of such global 
properties. Figure 5(b) illustrates a proposed manifestation of 
apparently distinct absorbers in the MHD wind discussed in this 
study. While it is beyond the scope of this work, we also note in 
a global viewpoint that the radial extent of magnetized winds is 
coupled to the disk magnetization (Murphy et al. 2010). 

In addition, there is an issue of the radial density dependence 
of these winds. As we discussed in the Introduction, the 
main motivating factor of our previous work, have been the 
observations of Behar (2009) that indicate Ah oc where s ~ 0 
and imply a LoS density profile n(r) oc r -( 2s+1 V0+i) ; 0 < s < 1 
with most values closer to s — 0 (the parameter s is related 
to that q of Equation (1) by 2q = (s + 2 )/(s + 1)). As we 
noted, this scaling suggests correlations between the column 
Ah, the velocity v and the ionization parameter £ of specific 
transitions at a given object (for instance Ah oc u 2l / (s+i \ v oc 
|(j+i)/ 2 ) More recently, observations, mainly of UFOs but 
also certain WA transitions, in an ensemble of (rather than 
a single) AGNs (Tombesi et al. 2010b, 2012a, T12c), have 
shown on the aggregate a preference for the value s — 1, 
implying a LoS density profile closer to n(r) oc r -3 '' 2 . At 
this point it would be of interest to note that models that 
include a combined accretion-outflow approach to the problem 
(Ferreira 1997; Ferreira & Casse 2004) are more restricted in 
their choice of parameters than those presented above and do 
provide additional restrictions to the values of s. So, “cold” 
flows (Ferreira 1997) are restricted to density profiles very close 
to n(r) oc r~ 3/2 , while models with “entropy injection” near 
the Alfven surface (Ferreira & Casse 2004) can yield density 
profiles close to n(r) oc r~ 1 . At this point we will reiterate that 
X-ray spectroscopy observations could provide the distinction 
between these extreme cases thus establishing a probe of the 
accretion-outflow physics of these systems. 

Another issue with the solutions we present is their reliance 
on self-similarity. The utility of this technique is not only 
mathematical, i.e., in allowing one to obtain solutions of 
the MHD equations; it is of importance because it provides 
solutions that span a multiple decades in radius. Solutions of 
this type are demanded by the data, i.e., the large range in 
ionization parameter spanned by the observed X-ray transitions, 
as discussed in Section 1 . Whether the simple functional form 
of the model parameters given above extends over the entire 
domain of validity of these solutions is an issue that should be 
answered observationally. Deviations from a single power law 
would indicate the presence of additional physics that will have 
to be included in such calculations. 

By their nature, self-similar solutions cannot cover the region 
near the flow symmetry axis and they must also terminate at 
some distance from the origin. As we have noted above, the 
near-axis high latitude region is much less dense and much 
more highly ionized than the lower latitude ones, so from the 
interpreting observed point of view of absorption features which 
is the main thrust of the present paper, it is not of significance. It 


is certainly of significance mathematically and for impacting the 
physics of jets, as the magnetic fields of this region are likely to 
thread the BH event horizon. Pressure balance requires that the 
pressure there match that at larger radii where the self-similar 
solutions are valid. This is possibly the pressure of relativistic 
particles that gives rise to the radio jets observed. Such models 
that combine a self-similar disk-wind solution with an MHD 
flow of limited lateral extent (i.e., one in which the magnetic 
field threads a star or a BH), have recently appeared in the 
literature (Matsakos et al. 2008, 2009; Porth & Fendt 2010). 
Again, it is very conceivable that the winds may consist of 
two components; an axial collimated jet of low matter density 
probably originating from the BH itself surrounded by an 
extended disk wind (e.g., Ferreira et al. 2010) also similar to 
our view depicted in Figure 5(b). It is of interest to note that 
the last of the above references finds flows with opening angles 
between 4° and 6°, which is not very different of those of models 
(A) and (D). 

At large radii, these solutions cannot be trusted beyond an 
equatorial distance larger than that of the underlying disk. This 
distance is not known. However, it cannot be larger than the BH 
influence radius, Rb ~ M^u/a 1 , i.e., the distance at which 
the disk Keplerian velocity matches the velocity dispersion 
<7 = (v 2 ) 1/2 of the AGN surrounding spheroid. However, given 
the relation between Mbh and a , namely (Mbh/ 10 8 M o ) ~ 
(er/200kms -1 ) 4 (Ferrarese & Merritt 2000), this distance is of 
the order of R B ~ 10 19 (M BH /10 8 ~ 10 19 M~ 1/2 cm. 

It is of interest to note that absorption feature observations, 
in conjunction with our models, could provide an independent 
estimate of this distance and therefore an independent measure 
of the M — a relation. 

Finally, one should bear in mind that density profiles n(r) oc 
r “ produce winds whose mass flux increases with distance for 
a A 3/2. Applying the profile with a — 1 and m = 0.1 used 
in our models to an accreting BH of M = 10 s M 0 , we obtain a 
total mass flux of 
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a value much higher than that needed to power its bolometric 
luminosity of ~10 45 erg/s. The corresponding (integrated) 
kinetic power 9 is given by £o° t cdl) = M (locaI) u 2 ut a 
Assuming r m /r a ~ 100, this yields £ou t obdl> ~ 10 44 er g S_1 th a t 
is indeed consistent with the estimated UFO power in 3C 111 
(Tombesi et al. 2010b, 2012a, T12c) providing a large impact 
on AGN feedback process at large scales. On the other hand, 
the cumulative power of UV /soft X-ray WAs may be able to 
reach this level as well (e.g., Crenshaw & Kraemer 2012). It 
is interesting to note that outflow rates, estimated by X-ray 
spectroscopy in BH binaries (Neilsen & Lee 2009; Neilsen et al. 
2012), were found to be much larger than the accretion rates 
needed to power their X-ray luminosity. 

On a closing note, although in this work we primarily focused 
on physical conditions and a global geometrical structure of 


9 Also see K12 for the discussion of radial dependence of the wind quantities. 
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MHD winds, we are currently simulating the corresponding 
spectral signatures of the absorption lines within this model (K. 
Fukumura et al. 2013, in preparation). To constrain individual 
wind parameters (such as F 0 , II and their dependency on m 
and 9), however, one may require numerical calculations of X- 
ray absorption feature with sets of MHD parameters and perform 
a / 2 -test in comparison with data; i.e., observationally constrain 
the wind parameters that would otherwise be inaccessible. While 
it is still challenging to constrain all the wind variables, one can 
still hope to restrict the range of wind parameters allowed by 
spectrum analysis. The capabilities of the upcoming mission 
Astro-H should be instrumental in performing these insightful 
observations. 

K.F. thanks an anonymous referee for improving the 
manuscript and M. Nakamura for his inspiring comments on 
the draft manuscript. E.B. is supported by grants from the Israel 
Science Foundation and from Israel’s Ministry of Science and 
Technology. 


APPENDIX 

A BRIEF COMPARISON OF OUR WIND 
PARAMETERS WITH BP82 WINDS 


We have computed and examined four magnetized wind solu- 
tions to steady-state, axisymmetric ideal MHD equations under 
Newtonian gravity (as in CL94, FKCB 10a, FKCB 10b, and K12) 
in the context of its potential role of outflow physics. Generic 
properties of MHD outflows considered in this study include 
mass loading and angular momentum as well as the global den- 
sity profile. We suggest that smaller angular momentum tends 
to help collimate the outflow poloidal structure while a lower 
plasma-to-magnetic flux ratio generally leads to a more efficient 
wind acceleration by the action of a large-scale magnetic fields 
(at least predominantly relative to other mechanisms such as 
radiation pressure). 

Here, we would like to show a viability of the sets of our wind 
parameters by recalling some of the fundamental definitions of 
the conserve quantities. As discussed in CL94, one can express 
F a in terms of some “known” quantities from Equation (7) 
combined with 

F(‘P) = 4“ 3/2 F 0 — , (Al) 

Vo 


to obtain 


rp . v p v 0 
F 0 = 4?rp — — 

Hp t ) 0 

~ 10 n (0.01 x (0.7c) 2 )/(10 4 ) 2 


O(10~ z ). (A2) 


assuming a typical set of values (for AGNs) of n a ~ 10 u) cm 3 
(Crenshaw et al. 2003 and references therein), B p ~ B„ ~ 
10 4 G, and v p ~ 0.01 v 0 at the innermost launching radius. This 
value is in consistence with our chosen values. For the rest of 
the parameters, we (this paper and CL94) are also convinced of 
consistency with our parameter sets in comparison with BP82 
winds as also discussed in CL94; i.e., one can make a one-to- 
one mapping between BP82 wind parameters and ours by noting 
that 


A. ** -H 0 /F 0 , (A3) 

where BP82’s characteristic value of X = 30 leads to —H a ~ 
0(1). From our boundary condition at the disk surface expressed 


in Equation (9) (see also CL94) one finds 

£2 C = 1 ~(F 0 + H 0 )V Z , o~0(l), (A4) 

where we have used F a and Fl a above together with V z o ~ 0.01 . 
Therefore, we are confident in our parameter choice. 
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